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ABSTRACT 

SN 2006tf is the third most luminous supernova (SN) discovered so far, after SN 2005ap and 
SN 2006gy. SN 2006tf is valuable because it provides a link between two regimes: (1) luminous type 
Iln supernovae powered by emission directly from interaction with circumstellar material (GSM), and 
(2) the most extremely luminous SNe where the GSM interaction is so optically thick that energy must 
diffuse out from an opaque shocked shell. As SN 2006tf evolves, it slowly transitions from the second 
to the first regime as the clumpy shell becomes more porous. This link suggests that the range in 
properties of the most luminous SNe is largely determined by the density and speed of hydrogen-rich 
material ejected shortly before they explode. The total energy radiated by SN 2006tf was at least 
7 X 10^° ergs. If the bulk of this luminosity came from the thermalization of shock kinetic energy, 
then the star needs to have ejected ~18 M© in the 4-8 yr before core collapse, and another 2-6 Mq 
in the decades before that. A Type la explosion is therefore excluded. From the Ha emission-line 
profile, we derive a blast-wave speed of 2,000 km s~^ that does not decelerate, and from the narrow P 
Cygni absorption from pre-shock gas we deduce that the progenitor's wind speed was ~190 km s~^. 
This is reminiscent of the wind speeds of luminous blue variables (LBVs) , but not of red supergiants 
or Wolf-Rayet stars. We propose that like SN 2006gy, SN 2006tf marked the death of a very massive 
star that retained a hydrogen envelope until the end of its life, and suffered extreme LBV-like mass 
loss in the decades before it exploded. 

Subject headings: circumstellar matter — stars: evolution — stars: mass loss — stars: winds, outflows 
— supernovae: individual (SN 2006tf) 



1. INTRODUCTION 

Massive stars that die as core-collapse supernovae 
(SNe) of Types II and Ib/c span a wide range in lumi- 
nosity, but they are typically fainter than the standard 
thermonuclear Type la events that mark the deaths of 
lower-mass stars. More-luminous counterexamples are 
usually observed as SNe of Type Iln, exhibiting rela- 
tively narrow [^ 1000 km s~^) emission lines of H in 
their spectra (Schlegel 1990; for a review of SN classifi- 
cation, see Filippenko 1997). These lines are generally 
attributed to shock interaction with dense circumstellar 
matter (GSM) rather than photospheric, high-velocity 
emission from the SN ejecta. They can show multiple 
components: a very narrow feature having typical widths 
up to a few hundred km s~^ attributed to emission from 
pre-shock GSM shed by the progenitor, as well as broader 
components with widths of a few thousand km s~^ aris- 
ing from dense post-shock gas. We will refer to these 
two components as "narrow" and "intermediate-width" 
lines, respectively. By comparison, broad lines in SN 
ejecta typically have widths of 10,000-20,000 km s^^. 

Sufficiently dense GSM can decelerate the blast wave 
and convert its bulk kinetic energy into X-rays and 
then visual radiation, thereby substantially increasing 
the bolometric luminosity of the SN. This is seen in lumi- 
nous SNe Iln and II-L having very strong Ha emission, 
such as SNe 1979G, 1986J, 1988Z, 1994W, 1997cy, and 



others (Branch et al. 1981; Filippenko 1991, 1997; Lei- 
bundgut et al. 1991; Turatto et al. 1993; SoUerman et al. 
1998; Chugai et al. 2004; Germany et al. 2000; Turatto 
et al. 2000; Benetti et al. 1998), as well as the "hybrid" 
Type la/IIn objects SN 2002ic (Hamuy et al. 2003; Wang 
et al. 2004; Wood-Vasey et al. 2004; Kotak et al. 2004; 
Benetti et al. 2006) and SN 2005gj (Prieto et al. 2007; 
Aldering et al. 2006). SNe 1997cy and 1999E may also be 
hybrid SNe of this type (Germany et al. 2000; Filippenko 
2000; Rigon et al. 2003; Hamuy et al. 2003; Wood-Vasey 
et al. 2004). 

Because it is the dense GSM that absorbs momentum 
and drains kinetic energy from the blast wave, more- 
luminous SNe Iln require progenitors with higher mass- 
loss rates or slower wind speeds. To account for some 
of the more-luminous SNe Iln, inferred progenitor mass- 
loss rates need to be of order 0.1 Mq yr^^ or higher 
(e.g., Chugai et al. 2004; Chugai & Danziger 2003; Smith 
2008). These extreme requirements point to episodic 
mass ejection reminiscent of the eruptions seen in -q Cari- 
nae and other luminous blue variables (LBVs), which 
have mass-loss rates of 0.01-1 Mq yr~^, far exceeding 
the limiting mass loss in line-driven stellar winds (Smith 
& Owocki 2006; Owocki et al. 2004). In SN Iln progen- 
itors, the reason for this extreme mass loss is not yet 
understood, but it must occur in the decades preceding 
core collapse because the blast wave can only reach radii 
out to a few 10^ AU during early phases of these SNe. 
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Fig. 1. — The top panel shows our BVRI light curves of SN 2006tf obtained with KAIT (see Table 1), and the bottom panel gives the 
apparent B — V color. The dotted line illustrates a least-squares fit to the color change with a slope of 0.0025 ±0.0003 mag d~^. Days 
since discovery are plotted because the explosion date is not known. 



Two recent events have pushed the limits on the re- 
quired physical conditions, challenging the standard pic- 
ture of CSM interaction: SN 2006gy (Ofek et al. 2007; 
Smith et al. 2007; Woosley et al. 2007; Smith & Mc- 
Cray 2007) radiated more energy in visual light than any 
other known SN, and SN 2005ap (Quimby et al. 2007a) 
appears to have had the brightest peak absolute magni- 
tude yet observed. To account for their extreme lumi- 
nosities with CSM interaction, the required progenitor 
mass- loss rates are of order 1 Mq yr~^. The problem is 
that in SN 2006gy, the weak Ha and X-ray emission im- 
ply mass-loss rates a factor of 10^-10'* smaller (Smith et 
al. 2007), while SN 2005ap shows no spectral signature of 
CSM interaction at all (i.e., it was not a SN Iln; Quimby 
et al. 2007a). One possible solution is that these SNe are 
indeed ultimately powered by CSM interaction, but this 
interaction occurs early in a highly opaque shell, forcing 
photons to diffuse out after the shock has already passed 
through it (Smith & McCray 2007) . An alternative view 
(Smith et al. 2007; Quimby et al. 2007a) is that these 
two SNe could be powered by radioactive decay from a 
large mass of ^^Ni synthesized in a pair-instability SN 
event (Barkat et al. 1967; Rakavy & Shaviv 1967; Bond, 
Arnett, & Carr 1984), but that hypothesis remains diffi- 
cult to prove or rule out (see Smith et al. 2008b), partly 
because the decay luminosity would also need to diffuse 
out through a massive envelope. 

In this paper we investigate yet another extremely lu- 
minous SN that is nearly as luminous as SNe 2006gy and 



2005ap and has some similarity to them, but which ex- 
hibits signs of stronger CSM interaction in the spectrum. 

SN 2006tf was discovered on 2006 Dec. 12 (UT dates 
are used throughout this paper) in the course of the 
Texas Supernova Search (Quimby et al. 2007b). In spec- 
tra obtained ^10 d after discovery it showed Type Iln 
signatures. SN 2006tf is located 0'.'2 east and 0'.'7 north 
of the center of an anonymous and faint galaxy with an 
apparent R magnitude of 20.68 (Quimby et al. 2007b), 
for which the metallicity is not known. Its redshift, 
z =- 0.074, indicates a distance of 308 Mpc, assuming 
Ho — 72 km s~^ Mpc~^. At that distance, the discov- 
ery magnitude of about rn^ — 16.7, which was roughly 
constant for the first 2 weeks after discovery, makes 
the peak absolute magnitude about —20.7, or roughly 
1.6x10^° L0 without correcting for extinction. 

This extremely high luminosity makes SN 2006tf one of 
the most luminous SNe ever discovered. Among known 
SNe, it is surpassed only by SN 2006gy and SN 2005ap. It 
is more luminous than very bright or long-lasting SNe Iln 
with strong CSM interaction like SNe 1979C, 1988Z, 
1994W, 1998S, and others, as well as the Type la/IIn 
"hybrid" objects mentioned earlier. Given its Type Iln 
spectrum, the extreme luminosity of SN 2006tf is likely 
to be powered by CSM interaction. With the lack of 
any Type Iln signatures in the spectrum of SN 2005ap 
(Quimby et al. 2007a), and the difficulty in reconciling 
the energy budget of SN 2006gy with its relatively weak 
CSM interaction signatures (Smith et al. 2007, 2008b), 



SN 2006tf Gets the Bronze 



SN 2006tf 



day 32 (j^y 54) 



7800 K 




4000 



5000 6000 7000 

Rest Wavelength (Angstroms) 



8000 



Fig. 2. — The top panel shows visual-wavelength spectra of SN 2006tf on days 32, 41, 66, and 95 since discovery, normalized to the red 
continuum level, and with constant offsets of +3 (d32), +2 (d41), and +1 (d66). The day 64 spectrum (obtained for spectropolarimetry) is 
plotted in red several times for comparison with the other epochs. The dotted blue curves show representative black bodies for comparison 
at the temperatures indicated. At the bottom, the day 64 spectrum of SN 2006tf (black) is compared to the day 93 spectrum of SN 2006gy 
(green), dereddened assuming E{B — V) = 0.723 mag (Smith et al. 2007), and to the day 49 spectrum of SN 1994W (blue) from our 
spectral database. All spectra of SN 2006tf have been corrected for E{B — V) = 0.027 mag. 

(KAIT; Filippenko et al. 2001) at Lick Observatory. Flat 
fielding and bias subtraction were processed automati- 
cally. Galaxy subtraction and differential photometry 
were done using the KAIT pipeline (Ganeshalingam et 
al., in prep.). Template images of the galaxy were ob- 
tained on 2008 Jan. 18 with KAIT, more than 1 yr af- 
ter discovery. Our late-time spectrum (see below) shows 
strong Ha emission and some faint continuum detected 
from the SN. Hence, there is some concern that our tem- 
plate image obtained on day 404 is partly contaminated 
by SN emission, rather than pure background galaxy 
emission. At this late time, we measure the magnitude 
at the position of the SN to be 20.1 in the R band. Even 
if this late-time magnitude is dominated by SN light, 
the consequent oversubtraction would have caused us to 
overestimate the true R magnitude by only 0.22 mag 
after day 150, which is not much larger than the photo- 
metric uncertainty at that time. 

Calibrations for the field were obtained on 8 photomet- 



the strong GSM signatures and comparably high lumi- 
nosity of SN 2006tf are of significant interest. Therefore, 
in this paper we consider the photometric and spectro- 
scopic evolution of SN 2006tf in detail. 

We present our optical observations in §2. The light 
curve and bolometric luminosity are discussed in §3, and 
in §4 we describe the general spectral evolution. A de- 
tailed look at the Ha and H/3 emission-line profiles and 
their evolution is given in §5, including a late-time spec- 
trum taken on day 445. Analysis of our spectropolarime- 
try of SN 2006tf is presented in §6. In §7 we estimate the 
most basic physical properties of SN 2006tf. Section 8 
summarizes our main results, discussed in the context of 
an integrated picture of SN 2006tf and other SNe pow- 
ered by GSM interaction. 

2. OBSERVATIONS 

We obtained optical photometry of SN 2006tf in BVRI 
using the 0.76 m Katzman Automatic Imaging Telescope 



Smith et al. 



TABLE 1 
KAIT Magnitudes of SN 2006tf=' 



MJD 


B 


err a 


V 


err^ 


R 


err^ 


/ 


err/ 


2454115.02 


17.739 


0.030 


17.199 


0.017 


16.923 


0.014 


16.707 


0.028 


2454118.99 


17.784 


0.036 


17.322 


0.028 


17.040 


0.021 


16.711 


0.052 


2454121.99 


17.852 


0.041 


17.300 


0.030 


17.032 


0.030 


16.767 


0.039 


2454125.00 


17.794 


0.025 


17.338 


0.019 


17.036 


0.014 


16.754 


0.031 


2454131.00 


17.773 


0.039 


17.379 


0.031 


17.123 


0.023 


16.720 


0.040 


2454133.94 






17.408 


0.079 






16.769 


0.129 


2454144.04 


18.183 


0.075 


17.601 


0.153 


17.177 


0.145 


16.730 


0.107 


2454146.97 


18.233 


0.075 






17.245 


0.072 


16.895 


0.035 


2454149.94 


18.184 


0.033 


17.517 


0.023 


17.179 


0.018 


16.892 


0.039 


2454161.98 


18.275 


0.370 


17.549 


0.080 


17.291 


0.054 


16.886 


0.070 


2454170.92 


18.461 


0.045 


17.837 


0.032 


17.455 


0.027 


16.987 


0.055 


2454175.90 


18.486 


0.047 


17.796 


0.031 


17.488 


0.028 


17.125 


0.045 


2454178.88 


18.598 


0.049 


17.858 


0.047 


17.571 


0.025 


17.119 


0.044 


2454181.86 


18.675 


0.058 


17.926 


0.028 


17.638 


0.031 


17.226 


0.037 


2454184.88 


18.627 


0.083 


18.115 


0.127 


17.614 


0.068 


17.235 


0.100 


2454188.86 


18.527 


0.092 


18.153 


0.064 


17.621 


0.037 


17.196 


0.059 


2454198.84 


18.814 


0.168 


18.231 


0.126 


17.849 


0.152 


17.336 


0.116 


2454201.80 


19.004 


0.086 


18.146 


0.045 


17.865 


0.029 


17.426 


0.041 


2454204.85 






18.252 


0.148 


17.889 


0.257 


17.549 


0.136 


2454207.82 


19.120 


0.073 


18.368 


0.062 


17.977 


0.043 


17.628 


0.058 


2454211.81 


19.186 


0.088 


18.348 


0.053 


18.012 


0.030 


17.647 


0.055 


2454225.81 


19.003 


0.097 


18.547 


0.064 


18.113 


0.046 


17.765 


0.087 


2454230.77 


19.281 


0.061 


18.626 


0.055 


18.239 


0.041 


17.867 


0.051 


2454235.79 


19.280 


0.120 


18.610 


0.122 


18.165 


0.080 


17.976 


0.071 


2454242.75 


19.710 


0.071 


18.798 


0.075 


18.401 


0.062 


17.967 


0.121 


2454249.71 


19.475 


0.137 


18.767 


0.137 


18.589 


0.160 


18.004 


0.192 


2454254.77 


19.442 


0.248 


18.767 


0.127 


18.374 


0.101 







Photometric uncertainties arc la. 



TABLE 2 
Spectroscopic Observations op SN 2006tf 



Date 


Day'' 


Instrument 


Range'' 
A 


A/AA 


EW(Ha) 
A 


F(Ha) 
10"''' erg s-l cm^2 


Ha/H/3'= 


N/T^ 
% 


2007 Jan. 13 


32 


LRIS/Keck I 


3500-8600 


1000 


-44(6) 


1.33(0.1) 


2.85 


17 


2007 Jan. 22 


41 


DEIMOS/Keck II 


4200-8900 


2000 


-66(5) 


1.81(0.17) 


3.67 


11.2 


2007 Feb. 14 


64 


LRIS/Keck I (pol) 


3500-8600 


700 


-72(5) 


2.19(0.16) 


3.48 


9.0 


2007 Feb. 16 


66 


DEIMOS/Keck II 


4300-6750 


8000 


-82(7) 


1.90(0.2) 


3.79 


8.1 


2007 Mar. 17 


95 


DEIMOS/Keck II 


4200-8900 


2000 


-194(11) 


3.39(0.24) 


4.84 


3.4 


2008 Mar. 1 


445 


ESI/Keck II 


4400-8800 


6000 


-1035(90) 


0.69(0.34) 


11.6 


3.6 



^ Days since discovery. The additional number of days since explosion could be as much as 50—70 d for SN 2006tf, if the behavior of SN 2006gy is 
relevant. Spectral range in rest wavelength.^ The Ha/H/3 flux ratio measured after the spectra are dereddened. Percentage contribution of the 
narrow component to the total Ha line flux. 



ric nights using both KAIT and the 1 m Nickel telescope 
at Lick Observatory. The uncertainty in our subtraction 
and photometry pipeline is estimated by injecting arti- 
ficial stars with the same magnitude and point-spread 
function of the SN into regions of comparable complex- 
ity in the original KAIT images and recovering them. 
The final uncertainty is taken to be the scatter in re- 
covering 20 artificial stars added in quadrature with the 
calibration error. 

We take day zero to be MJD = 2,454,081.98, the dis- 
covery date 2006 Dec. 12 (Quimby et al. 2007b). Since 
this is not the explosion date, there will be a source of 
continual ambiguity when comparing SN 2006tf to other 
SNe. The explosion date is not known, but if the light- 
curve evolution was similar to that of SN 2006gy, the 
explosion of SN 2006tf could conceivably have been long 
before the time of discovery. Figure [T] shows the KAIT 
magnitudes, which are listed in Table 1. 

During the main peak of the light curve (within 200 d 



of discovery), we obtained visual- wavelength spectra of 
SN 2006tf on five separate dates at the Keck Observatory 
using the low resolution imaging spectrometer (LRIS; 
Oke et al. 1995) and the Deep Imaging Multi-Object 
Spectrograph (DEIMOS; Faber et al. 2003), with the ob- 
servations summarized in Table 2. The epochs for these 
four spectra are days 32, 41, 64, 66, and 95 after dis- 
covery. The spectra were all obtained with the long slit 
oriented at the parallactic angle (Filippenko 1982), and 
were reduced using standard techniques. Wavelengths 
were corrected for redshift z = 0.074 so that the nar- 
row Ha emission line was at the proper rest wavelength. 
Given the very weak Na I D absorption and blue color, we 
made no correction for reddening and extinction beyond 
the relatively small Galactic values of E{B — V) = 0.027 
mag and Ar = 0.062 mag (Schlegel et al. 1998). The 
resulting spectra are shown in Figure [2l 

We also obtained one late-time optical spectrum on 
2008 March 1, day 445 after discovery, using the Echelle 



SN 2006tf Gets the Bronze 



Spectrograph and Imager (ESI; Sheinis et al. 2002) on 
the Keck II telescope (see Table 2). A 1200 s exposure 
was obtained in echellete mode, with the O'.'75-wide slit 
oriented at position angle 45° . Because of relatively poor 
seeing of ~1'.'5 (wider than the slit) and possible light 
cirrus during the night, the flux calibration is uncertain, 
and we conservatively adopt a factor of 2 uncertainty in 
absolute flux (i.e.. Ha in Table 2). This is larger than the 
line-measurement uncertainty, and is hard to quantify. 
The late-time spectrum is discussed in §5.4. 

Thus, we have two independent estimates of the late- 
time R magnitude. As we mentioned earlier, our KAIT 
template image shows a magnitude of 20.1, or L = 
7.5 X 10^ Lq, at the position of the SN. The true SN 
light could be substantially fainter than this if back- 
ground galaxy light is important, but it cannot be much 
brighter (note that Ha also contaminates this measure- 
ment). From our late-time spectrum on day 445, in 
which we subtracted nearby background and still de- 
tected faint continuum emission from the SN, we measure 
Fx = 1.4 ergs s""'^ cm~^ A^^ for the continuum level at 
red wavelengths, or an apparent R magnitude of 20.3. 
This corresponds to a late-time bolometric luminosity 
oi L — 6.3 X 10^ L0. Even though our uncertainty in 
the flux is a factor of two because of sky conditions, the 
agreement with the late-time KAIT magnitude to better 
precision than that is reassuring. 

Of the five epochs of spectra during the first 100 d, 
we obtained one epoch of spectropolarimetric observa- 
tions of SN 2006tf on 2007 Feb. 14.54 (day 64) using the 
polarimetcr unit^ of LRIS. The object was observed for 
1000 s in each of the four rotation angles of the half- 
wave plate retarder. The total observed spectral range 
was 3350-9240 A, and the l'.'5-wide slit gave a spectral 
resolution of ^^5 A in the blue and ~9 A in the red. 
The spectropolarimetric reductions followed the proce- 
dure outlined by Miller et al. (1988) and implemented 
by Leonard et al. (2001). Observations of the polar- 
ized standard stars HD 19820 and HD 161056 (Schmidt 
et al. 1992) gave consistent results for the zero point 
of the instrument. Instrumental polarization was neg- 
ligible as both of the observed null polarization stan- 
dards showed <0.1% polarization. Two foreground stars 
(BD -hll 2491, BD -hl2 2510), both within 20' of the SN, 
were observed to constrain any potential Galactic compo- 
nent of interstellar polarization (ISP) due to dust grains 
along the line of sight. Both stars showed less than 0.1% 
polarization, indicating that the Galactic component of 
ISP is also negligible, consistent with the low dust col- 
umn of E{B -V)^ 0.027 mag at the high Galactic 
latitude of SN 2006tf (6 =74°). 

3. THE LIGHT CURVE AND BOLOMETRIC LUMINOSITY 

Our BVRI photometry of SN 2006tf obtained with 
KAIT is presented in Figure [H Two features are im- 
mediately apparent. The first is the remarkably slow de- 
cline in the light curve, with roughly 0.01 mag d~^ in the 
R band, reminiscent of the slow evolution in SN 2006gy 
as well as the very slow decline in strongly interacting 
SNe Iln like SNe 1988Z, 1997cy, and 2005gj. The decline 
rate also matches the expected decay rate of ^^Co, and 

^ See http : //alamoana . keck . hawaii . edu/inst/lr is/polarimeter / 1 
for the online polarimeter manual by Cohen (2UU5J. 
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Fig. 3. — The absolute _R-band light curve of SN 2006tf compared 
to that of several other luminous or strongly interacting SNe II 
and Iln. SN 2005ap is a Type II (not a Iln) from Quimby et 
al. {2007a) and SN 1999em is a normal Type II-P SN (Leonard 
et al. 2002). All others shown are Type Iln: SN 2006gy data 
are from Smith et al. (2007), SN 1994W data from Sollerman et 
al. (1998), SN 1988Z from Turatto et al. (1993), and SN 1998S 
from Fassia et al. (2000). SN 1997cy (Germany et al. 2000) is an 
extremely luminous SN Iln possibly associated with GRB 970514, 
or alternatively, a possible "hybrid" Type la/IIn like SNe 2002ic 
(Wood-Vasey et al. 2004; Hamuy et al. 2003) and 2005gj (Prieto 
et al. 2007) that are still difficult to understand (see also Benetti 
et al. 2006). When different values of Hq were used by the above 
authors, the absolute magnitudes have been adjusted to correspond 
to our adopted value of 72 km s~^ Mpc~^. 

this is discussed in more detail below. The second feature 
is the gradual change in color. The lower panel in Fig- 
ure [T] shows the observed B — V color, which gets redder 
at a rate of ~ 0.0025 ± 0.0003 mag d"^ This slow evo- 
lution to redder colors is in agreement with our spectra 
in Figure [51 in which the continuum slope gradually and 
monotonically evolves to lower temperatures. The rela- 
tively weak color changes, with temperatures that imply 
little or no bolometric correction, mean that we can take 
the absolute R magnitude (corrected for small Galactic 
extinction) as a proxy for the bolometric magnitude, to 
compare SN 2006tf with other luminous SNe. 

Thus, in Figure [3] we compare the absolute R mag of 
SN 2006tf to that of several other notable SNe Il/IIn. 
Near its peak at the time of discovery, SN 2006tf 
was the third most luminous SN observed so far af- 
ter SN 2005ap and SN 2006gy. However, SNe 2005gj, 
2002ic, and 1997cy are not far behind. All of these are 
SNe II. Interestingly, though, all except the most lumi- 
nous, SN 2005ap, were SNe Iln with prominent and rel- 
atively narrow Ha indicating strong GSM interaction.^ 
The strength of the SN Iln signature may be an impor- 
tant clue to the nature of the most luminous SNe: As 
we noted in a previous paper (Smith ct al. 2007), while 
SN 2006gy was definitely a Type Iln with an Ha pro- 

2 The true natures of SNe 2002ic, 2005gj, and 1997cy are still 
controversial: SNe la exploding in dense H envelopes or perhaps 
the deaths of massive stars; see Benetti et al. (2006). 
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file that indicated some CSM interaction, the strength 
of that interaction indicated by the Ha and X-ray lumi- 
nosity seemed far too weak to be compatible with the 
power in continuum light. SN 2005ap was an even more 
extreme case, with no detectable sign of a narrow Ha 
feature generated by CSM interaction. SN 2006tf may 
be an important link, as its continuum luminosity and 
strong Ha are intermediate between these two extreme 
cases and other strongly interacting SNe. This will be a 
recurring theme in our investigation. 

SN 2006tf was apparently discovered at, or perhaps 
somewhat after, its time of peak luminosity. The ex- 
plosion date and the photometric behavior during the 
rise to maximum are not known. This should be kept 
in mind when interpreting Figure [3l because the trans- 
lation of SN 2006tf on the time axis is necessarily quite 
uncertain. Judging by its similarity to SN 2006gy, for 
example, SN 2006tf could potentially be shifted as much 
as 50-70 d to the right to be more correctly compared 
with other SNe plotted here. Indeed, its high luminosity 
requires a large emitting radius of order 300 AU. That, 
combined with its maximum observed expansion speeds 
of order 7500 km s~^ seen in the broad component of the 
Ha line (see below), argue that it may have had quite a 
long rise time of order 50-100 d. Judging by SN 2GG2ic 
and SN 2005gj, both of which were caught during the rise 
to maximum, SN 2006tf should be shifted to the right in 
Figure [3] by at least 20 d. 

In any case, SN 2006tf clearly has an extremely high 
peak luminosity and a very long duration, with no sign, 
as yet, of a sharp decline like that of SNe 1994W or 
2005ap (Fig. [3]). As such, it demands an explanation 
for its high fluence that far exceeds that of normal core- 
collapse SNe. In general, such longevity at high lumi- 
nosity implicates a large mass of emitting material, as 
we shall see below. Although less extreme, the energy 
demands of SN 2006tf rival those of SN 2006gy, and 
considerations for the possible energy source are simi- 
lar. More in-depth considerations of the possible energy 
sources and associated difficulties can be found in our 
previous papers on SN 2006gy (Smith et al. 2007, 2008b; 
Smith & McCray 2007). The two sources of luminos- 
ity we consider here are ^^Co decay from a large initial 
mass of '^^Ni, or shock-deposited kinetic energy diffusing 
out of an opaque envelope after being thermalized and 
converted to visual light. 

3.1. Radioactive Decay? 

One possible source for a long lasting, extremely high 
SN luminosity is radioactive decay from a large initial 
mass of ^^Ni. This might be synthesized in a "hyper- 
nova" or coUapsar as has been proposed in the case of 
SN 1997cy where 2.3-2.6 Mq of ^^Ni was derived (Ger- 
many et al. 2000; Turatto et al. 2000); in that case, it is 
interesting that the light curve of SN 1997cy is quite sim- 
ilar to that of SN 2006tf. Alternatively, a large mass of 
^^Ni could potentially be synthesized in a pair-instability 
SN explosion (Bond et al. 1984; Barkat et al. 1967; 
Rakavy & Shaviv 1967), as had been noted as a pos- 
sibility for SN 2006gy (Smith et al. 2007) and SN 2005ap 
(Quimby et al. 2007a). 

Figure |3] shows the bolometric luminosity light curve 
of SN 2006tf, derived from the absolute R magnitude 
(adopting An = 0.062 mag). The thick gray line shows 
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Fig. 4.— The absolute R magnitudes of SN 2006tf from Figure [3] 
are converted to luminosities in Lq; we also include the uncer- 
tain late-time magnitudes mentioned in §2. These are compared 
to a simple photon diffusion model (Smith & McCray 2007) and 
the decline rate of ^^Co decay. The diffusion model (long dashes) 
gives an acceptable fit to the light curve for a diffusion timescale 
*diff = 160 d and a CSM shell of ~18 Mq, while radioactive decay 
(thiols gray line) would require ^4.5 Mq of ^^Ni to power the same 
luminosity (see text). We also show a CSM interaction model for 
a progenitor star of 110 Mq and a circumstellar envelope mass of 
~25 Mq (Woosley et al. 2007). 



the ^^Co luminosity and decay rate one would expect 
from an initial ^^Ni mass of about 4.5 M©.^ Although it 
does not match the light curve precisely, radioactive de- 
cay does reproduce the average decline rate of SN 2006tf 
quite well. It gives as good an account of the light 
curve as the 2.3-2.6 Mq of ^^Ni proposed in the case 
of SN 1997cy (Germany et al. 2000; Turatto et al. 2000). 

Furthermore, radioactive decay makes a clear predic- 
tion for the late-time luminosity: it should follow the 
^^Co decay rate. As we noted earlier, on day 404, a 
source at the position of SN 2006tf had an R magnitude 
of roughly 20.1 (or somewhat fainter) in the template im- 
age that we used to subtract background light and cal- 
ibrate our KAIT photometry in Figure [TJ Also, from 
the continuum level in a late-time spectrum obtained 
on day 445, we estimate an apparent R magnitude of 
20.3±0.75, as noted above. The late-time luminosity pre- 
dicted by the same ^^Co decay from the earlier light curve 
is roughly consistent with these late-time measurements. 
The decay model actually underpredicts the luminosity, 
but so do the CSM interaction models in Figure [H and 
hence any model requires us to invoke some additional 
CSM interaction at late times. 

Thus, we cannot easily discount ^^Co decay as a po- 
tentially significant source of luminosity for SN 2006tf, 
especially if some portion of the luminosity is supplied 
by CSM interaction as indicated by the spectral char- 
acteristics. This was also the case for SN 1997cy and 
SN 2006gy. Although we cannot confidently rule out ra- 
dioactive decay, it seems more likely to us that SN 2006tf 
is dominated by extreme CSM interaction. 

•^ The necessary ^®Ni mass is very uncertain. It could be 7 Mq 
if our adopted i = is actually ^SO d after the explosion date, 
or it could be arbitrarily less if CSM interaction contributes some 
substantial fraction of the luminosity. 
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3.2. Photon Diffusion From a Massive Shocked 
Envelope ? 

A better match to the main light-curve peak of 
SN 2006tf is attained with an opaque sheU-shocked model 
like that which Smith & McCray (2007) proposed for the 
main light-curve peak of SN 2006gy (Fig. |4]). In such a 
model, the extreme luminosity is a product of the large 
initial radius of an optically thick CSM ejecta envelope, 
throughout which shock kinetic energy has been thermal- 
ized efficiently If the envelope is highly opaque, radia- 
tion must diffuse out rather than being emitted directly 
by the post-shock cooling zone. This is true no matter 
what the ultimate source of energy deposition is, whether 
it be shock energy or radioactive decay, or both. 

The main ingredients that determine the shape of the 
diffusion light curve are the characteristic time for pho- 
ton diffusion, tdiff = tR/c, and the expansion time, texp 
(see Smith & McCray 2007; Falk & Arnett 1973, 1977). 
These can be expressed as idiff ~ 23(A//i?i5) d, with M 
being the envelope mass in M© and i?i5 being the char- 
acteristic initial shell radius in units of 10^^ cm. The ex- 
pansion timescale is just ioxp ~ (Ai?/ucxp) for constant 
expansion speed. Our model (dashed curve in Fig. 3]) is a 
simple analytic model with tdis — 160 d, i?i5 — 2.7, and 
M=18 M0. For comparison. Figure H] also shows a one- 
dimensional (1-D) numerical model of a similar physical 
situation from Woosley et al. (2007), with a 25 Mq shell 
ejected by a 110 M0 progenitor star about 5 yr before the 
SN. This model was intended for SN 2006gy, but seems 
even more applicable to SN 2006tf.^ 

The diffusion curve in Figure |4] provides an excellent 
fit to the observed data by adopting tdiff — icxp — 160 d 
(maximum luminosity for a given mass is achieved when 
these two are equal, minimizing adiabatic losses) . If the 
characteristic expansion speed of the blast wave is about 
2,000 km s'-^, as indicated by the intermediate- width Ha 
emission (see below), then i?i5 = 2.7 or about 180 AU 
for toxp — 160 d. This is the initial radius of the shocked 
ejecta envelope, which will grow with time. Thus, it is 
in approximate agreement with the somewhat larger ra- 
dius of ~300 AU that we would nominally expect for 
a ~7800 K black body emitting the peak luminosity of 
SN 2006tf (see below). For iaiff = 160 d, the required 
envelope mass is ~18 Mq. These values are similar 
to, but even more extreme than, those of SN 2006gy: 
idiff = 70 d, i? = 160 AU, and M = 10 M©, respectively 
(Smith & McCray 2007). As in that case, this radius 
is far too large to be the actual hydrostatic radius of a 
red supergiant progenitor, and the CSM-wind expansion 
speed of 190 km s~^ noted earlier is much too fast for that 
interpretation as well. Instead, this large radius probably 
represents the pseudo-photosphere of an opaque pre-SN 
ejecta shell that is not bound to the star, having been 
ejected ~4 yr prior to the SN if it has been moving at 
'--^190 km s^^. Containing 18 M© within a radius of 180 
AU would give at least r « 400. 

In this model, a self-consistent explanation for the 
lower peak luminosity and slower decline rate of 
SN 2006tf as compared to SN 2006gy is that the for- 
ward shock of SN 2006tf has a slower expansion speed of 

** Woosley et al. (2007) note that the peaks in the light curve 
may be an artifact of their 1-D models, which they expect would 
be smoothed out in reality. 



2,000 km s-\ compared to 4,000 km s^^ for SN 2006gy 
(Smith et al. 2007). The slower expansion speed causes a 
lower peak luminosity and longer duration, because the 
same amount of thermal energy will take a longer time 
to leak out of the envelope and the radius is smaller at a 
given time after explosion. 

Although a simple model like this can fit the main 
light-curve peak, it is not necessarily a unique explana- 
tion, and both types of CSM interaction models (Smith 
& McCray 2007; Woosley et al. 2007) fall short of the ob- 
served late-time luminosity (Fig. [4]). As with SN 2006gy, 
there must be some additional contribution of emission 
directly from the ongoing, more optically thin CSM in- 
teraction region in order to explain the strong Ha emis- 
sion. This occurs subsequent to the shock wave passing 
through the massive opaque shell, before which time the 
post-shock Ha emission would not be seen. 

3.3. Total Radiated Energy 

Regardless of which interpretation is correct, the pho- 
ton diffusion model in Figure |4] gives an accurate phe- 
nomenological fit to the bolometric luminosity, so it can 
be used to measure the total radiated energy. Integrat- 
ing this curve in Figure [5] from day to 180, we find 
that the total energy radiated in visual light during the 
time of our photometric monitoring was (6.2 ±0.3) x 10^° 
ergs, if we assume no bolometric correction to the abso- 
lute R magnitude. If we include the late-time tail, this 
value rises to about 7.0 x 10^° ergs, and it could be in- 
creased even further if the time of discovery occurred 
significantly after the actual time of explosion. This is 
likely to be the case if SN 2006tf had a slow rise time 
comparable to SN 2006gy, or even a faster rise akin to 
SNe 2002ic and 2005gj, as we noted earlier. 

It is interesting to note that an 18 M0 shell moving at 
2000 km s^^ contains the same amount of kinetic energy 
of about 7 X 10^° ergs. Rough equipartition of thermal 
and kinetic energy is not surprising in this model (see 
Smith & McCray 2007). 

In any case, the total amount of energy radiated in 
visual light is almost 10^^ ergs, or about half that of 
SN 2006gy (Smith et al. 2007). This amount of radiated 
energy must drain the reservoir of total available kinetic 
energy. Since the SN showed no sign of deceleration, 
this suggests either (1) that the CSM interaction we see 
now is not the major power source for the light curve, 
occurring subsequent to the shock passing through the 
opaque envelope mentioned earlier, or (2) that SN 2006tf 
marked an unusually energetic explosion well in excess of 
the canonical kinetic energy of a SN, as was the case with 
SN 2006gy (see Smith et al. 2007 for further discussion) . 

4. THE VISUAL-WAVELENGTH SPECTRUM AND 
CONTINUUM 

The top panel in Figure [2] shows our visual- wavelength 
spectra of SN 2006tf on the four dates during the decline 
from the main peak of the light curve, normalized to the 
level of the red continuum fiux and offset by constant 
values as noted in the caption. The epoch on day 64 
is repeated in red for comparison to the other 4 epochs 
(this epoch has the highest signal-to-noise ratio because 
it was obtained for spectropolarimctry; see §6). These 
spectra have been corrected for Galactic reddening of 
E{B -V)^ 0.027 mag. 
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One of the most striking changes with time is in the 
relative strength of Ha, to be discussed in detail in the 
following section. Here we focus mainly on the contin- 
uum shape and overall properties of the spectra. 

In general, Figure [2] reinforces our earlier conclusion 
from photometry that SN 2006tf shows mild change in 
its color and continuum shape during the first 100 d 
after discovery, with correspondingly little evolution in 
the character of the spectrum. Matching the continuum 
with black bodies (dotted blue curves in the top panel of 
Fig. [2]), the characteristic temperature drops monoton- 
ically from 7800 K to 6300 K between days 32 and 95, 
consistent with the steady change in B — V (Fig. [1]) . 
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Fig. 5.— Observed line profiles of Ha and H/3 in SN 2006tf. Ha 
has been normalized with a fit to the observed continuum level. 
H/3 for each date is normalized and then scaled arbitrarily for com- 
parison with the Ha profile. Day 64 is not shown here or in the 
next few figures because it is only 2 d before day 66 and has lower 
spectral resolution. 

Interestingly, the slow drop in temperature and charac- 
teristic line-blanketing features shortward of 4500 A are 
atypical for strongly interacting SNe Iln, which tend to 
show very blue continua at higher temperatures that re- 
main nearly constant if the blast- wave speed is constant. 



Instead, the slow decline in temperatures around 7000 K 
and pronounced line blanketing are more reminiscent of 
normal SNe H-P, where photons are diffusing out of the 
expanding SN ejecta through which the recombination 
photosphere is receding. Since SN 2006tf is far too lu- 
minous for this emission to arise in normal SN ejecta, 
this is another argument favoring the shell-shocked pho- 
ton diffusion model that we mentioned earlier, and the 
sustained high luminosity means that this shell must be 
very massive. SN 2006tf is distinguished from SNe II-P, 
however, in that its observed constant expansion speed 
indicates that this process must occur in a thin expand- 
ing shell rather than the geometrically thick envelope 
with a large velocity gradient in SNe II-P. 

In the bottom panel of Figure [2] we compare the spec- 
trum of SN 2006tf (black) to that of two other well- 
studied SNe Iln from our spectral database. SN 1994W 
(blue) is a case where the spectrum and luminosity are 
dominated by ongoing GSM interaction with a contribu- 
tion from diffusion of radiation from the SN ejecta or a 
shocked opaque envelope (Chugai et al. 2004), whereas 
in SN 2006gy (green) , most of the continuum luminosity 
arises from photon diffusion from a previously shocked 
opaque shell (Smith & McCray 2007).^ The continuum 
shape of SN 2006gy matches that of SN 2006tf quite 
well, as noted by Smith et al. (2007), but the Ha line 
of SN 2006tf is much stronger relative to its continuum. 
Taking into account both the continuum shape and the 
strong Balmer-line emission, we find that a linear combi- 
nation of the spectra of both SN 1994W and SN 2006gy 
gives a suitable approximation to SN 2006tf. This sug- 
gests that ongoing CSM interaction, as compared to dif- 
fusion, contributes a larger fraction of the observed lu- 
minosity in SN 2006tf than in SN 2006gy. Furthermore, 
the increase in the relative strength of the broad compo- 
nent of Ha with time implies that CSM interaction con- 
tributes a larger fraction of the luminosity at late times 
as the continuum source fades. This makes a strong case 
that SN 2006tf is a valuable example of the transition re- 
gion between opaque shocked shells, like in SNe 2006gy 
and 2005ap, and progenitors with lower mass-loss rates 
where CSM-interaction is seen directly, like in SN 1988Z. 

Finally, we note that aside from its Ha emission, the 
spectrum of SN 2006tf does not resemble the "hybrid" 
SNe Iln/Ia objects such as SNe 2002ic and 2005gj at any 
epoch, even though their light-curve shapes and absolute 
magnitudes are similar at early times. We detect no ev- 
idence for features from an underlying SN la spectrum. 
Additionally, the required envelope mass of 18 Mq and 
total radiated energy of almost 10^^ ergs make a SN la- 
powered interpretation unlikely for SN 2006tf. 

5. EVOLUTION OF Ha AND H/3 EMISSION 

Figure[5]shows both Ha and H/3 at all four epochs, with 
the Ha flux normalized to the nearby continuum, and 
the H/3 profile scaled for comparison. It is evident that 
H/3 is more asymmetric than Ha, showing a persistent 
deficit in the blueshifted emission wing at all epochs after 
day 32, due to high optical depths and the stronger self- 

^ The fact that we are comparing a day 64 spectrum of SN 2006tf 
to a day 93 spectrum of SN 2006gy is not necessarily cause for 
concern, since one is the date since discovery, and the other is 
the much longer time since the presumed explosion date. They 
represent comparable times after peak luminosity. 
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Fig. 6. — Observed Ha profiles on the four dates indicated, 
normalized and superposed over one another. The first three dates 
are plotted without further scaling, but the last one on day 95 has 
been adjusted by dividing the normalized intensity by 3 and adding 
0.66 to compare the line shape with the other dates. Velocities at 
which [N ll] AA6548, 6583 would be observed are indicated. 

absorption in H/3. 

The intermediate-width component of Ha at roughly 
±2000 km s~^ is discussed in §5.1, while broad emission 
at relatively high speeds is discussed in §5.2. In §5.3 we 
describe the profiles of the narrowest component from the 
unshocked CSM, in §5.4 we present the late-time spec- 
trum, and in §5.5 we discuss the behavior of the line 
luminosity with time. 

5.1. Post-Shock Ha and HP Profiles 

Figure [S] shows the intermediate-width Ha profile on 
days 32, 41, 66, and 95. This component at ±2000 km 
s""'^ dominates the appearance of the spectrum and ac- 
counts for most of the flux in the Ha line. The Ha profile 
appears to be identical on days 32 and 41, and the red 
wing of the line maintains the same shape at all epochs. 
However, we see a stronger, broad, blucshifted wing in 
the day 66 and 95 spectra. This enhanced blueshifted 
emission wing makes the line appear asymmetric. How- 
ever, we suspect that it is not related directly to the 
emission from post-shock gas, but instead is due to an 
underlying broad component that is seen as the shell be- 
comes more effectively optically thin, as discussed in the 
next section. Ignoring this broad blueshifted wing for 
now, the intermediate-width component of Ha is sym- 
metric and centrally peaked in SN 2006tf, typical of the 
Ha profiles that define SNe Hn. 

The speed of the blast wave expanding into the CSM 
can be inferred most clearly from the width and shape 
of this intermediate-width Ha profile. Figure [7] shows 
the Ha profile of SN 2006tf on day 41, when it is most 
symmetric, and compares it to the observed broad Ha 
profile of SN 2006gy from Smith et al. (2007). The reason 
this is useful is because the broad P Cygni absorption 
feature in SN 2006gy has a sharp blue edge at —4,000 km 
s^^ that provides a good measure of the likely speed of 
the blast wave (Smith et al. 2007), whereas the profile of 
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Fig. 7. — Same as Figure[6]but comparing the day 41 spectrum 
of SN 2006tf (solid black) to the day 96 Ha profile of SN 2006gy 
(shaded gray, from Smith et al. 2007; the flux has been scaled 
upward because the Ha emission of SN 2006gy is much weaker 
relative to its continuum level). The profile of SN 2006gy is no- 
ticeably broader than that of SN 2006tf. If we expand the velocity 
scale of SN 2006tf by a factor of 2 (dashed line), then the red side 
of the broad Ha profile matches that of SN 2006gy quite well. This 
suggests that the blast-wave speed of SN 2006tf is roughly half as 
fast as that of SN 2006gy. 

SN 2006tf has no sharp blue edge to provide a convenient 
measure of its true speed. 

The narrower Ha profile of SN 2006tf is an obvious 
clue that its blast-wave speed is significantly slower than 
that of SN 2006gy. If we artificially stretch the velocity 
scale of the day 41 Ha profile of SN 2006tf by a factor 
of 2, shown by the dotted line in Figure [3 then the red 
wing of the line matches that of SN 2006gy almost per- 
fectly. (A comparison on the blue side does not work 
because of P Cygni absorption in SN 2006gy.) This pro- 
vides a strong case that the blast-wave speed in SN 2006tf 
is about half that of SN 2006gy. Therefore, we adopt 
2,000 km s~^ as the characteristic blast-wave speed for 
SN 2006tf, defining the speed at which the CSM is swept 
up by the forward shock. This also matches the full width 
at half-maximum intensity (FWHM) of the intermediate 
component of the Ha line in SN 2006tf. Judging from 
the constant red wing (Fig. [S]), this speed remains un- 
changed for at least the first 100 d, during a time when 
the SN is radiating away ~10^^ ergs in visual light. 

5.2. Underlying Broad Component 

As noted earlier, comparing the Ha profiles at all four 
epochs (Fig. [5]) shows that the red side of the line changes 
very little, but there is extra emission that develops with 
increasing strength on the broad blueshifted wing. One 
would normally attribute the apparent flux deficit in the 
red wing to high optical depth that blocks the far side 
of the SN, but this hypothesis predicts a trend opposite 
to the one seen; i.e., if the asymmetry were due to high 
optical depths alone, then one would normally expect the 
broad emission to become more symmetric with time. 

Figure [8] investigates the blueshifted wings of broad- 
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Fig. 8. — {Top panel) Same as Figure |6l but with an expanded 
velocity scale and smaller vertical range to show the behavior of 
the broad blueshifted wings of Ha. The lower tracings show the 
day 66 and 95 residual profiles after subtracting off the symmetric 
day 41 profile. Residual fast blueshifted emission is seen out to 
roughly —7500 km s~^. (Bottom panel) Line profiles of O I A7774 
on day 64 (green) and 95 (orange), and the feature at ^5900 A 
(colors for different epochs are the same as the top panel), which 
could potentially be either He I A5876 or Na I A5892. It is plotted 
here as if it were He I A5876 (we have subtracted 0.3 from the 
normalized flux for display purposes), because that is the correct 
identiflcation for its narrow P Cygni component. 

line profiles in more detail. The profile of this "extra" 
emission on days 66 and 95 is illustrated in Figure [5] by 
subtracting the symmetric day 41 profile from the latter 
two epochs. The residual emission shows a broad com- 
ponent to Ha at blueshifted velocities, with its fastest 
speeds reaching about —7500 km s~^, and little or no 
corresponding emission on the red side of the line. 

Two possible explanations for this fast blueshifted 
emission are that it is due to electron scattering or to 
rapidly moving material. If it were due primarily to elec- 
tron scattering, though, we should see a component of 
comparable strength on the red side of the line as well. 

In fact, the explanation for this fast blueshifted Ha 
emission is not likely to be electron scattering alone, be- 
cause we also see evidence for fast-moving blueshifted 
material in absorption in some other species. The bot- 
tom panel of Figure [5] shows broad profiles of O i A7774 
and He i A5876,^ both of which indicate weak absorp- 
tion features with ill-defined blue absorption edges out 
to about —7000 or —8000 km s^^, matching the speeds 
at which we see enhanced blueshifted Ha emission. We 

^ We have labeled this feature as He I A5876 because that places 
its narrow P Cygni feature at the correct velocity, but the underly- 
ing broad feature could also plausibly be identified as Na I A5895. 



also note that the relative strength of the broad emission 
component of He i A5876 increases with time. 

Since the speeds up to about 7500 km s~^ are well in 
excess of the blast-wave speed we derived in the previ- 
ous section, the broad emission probably arises in gas 
that has not yet reached the reverse shock of the CSM 
interaction region. This expanding material just interior 
to the reverse shock may be ionized by the backward- 
propagating radiation from the post-shock gas (Cheva- 
lier & Fransson 1994). The severe asymmetry of the line 
means that we are only seeing the near side of the ejecta 
because of dust formation, high optical depths in the 
ejecta, or global asymmetry. 

The fact that these broad features are also seen in 
blueshifted absorption in He I A5876 and O i A7774 pro- 
vides a critical clue to their origin and to the physical 
processes at work in this SN. To be moving this fast, the 
material must have not yet reached the reverse shock, 
but to be seen in absorption, it must also have a back- 
ground light source. Thus, even as late as day 95, some 
central source (i.e., the inner SN ejecta?) still contributes 
a non-negligible fraction of the continuum luminosity. 

What constitutes a "non-negligible" fraction in this 
case? The background continuum source must contribute 
at least 10% of the luminosity at red wavelengths, be- 
cause that is roughly the depth of the absorption in O i 
A7774. This requires the source in question to have a lu- 
minosity of roughly 8 x 10* Lq (or Mr « —17.4 mag) on 
day 95, which is about 1 mag more luminous than a nor- 
mal SN II-P at that same time and even somewhat more 
luminous than a "hypernova" like SN 1998bw at that 
epoch (Galama et al. 1998). Thus, we must speculate 
that either the underlying SN ejecta of SN 2006tf were 
unusually luminous independent of its strong CSM inter- 
action, or that perhaps the strong CSM interaction has 
somehow rejuvenated the underlying ejecta with inward- 
propagating radiation. The alternative is that we would 
need the 7500 km s"-^ to get out ahead of the 2000 km 
s^^ shell, which seems unlikely without radical depar- 
tures from spherical symmetry. 

In these broad features, we have another clue that 
SN 2006tf is a more luminous analog of SN 1988Z, which 
also showed underlying broad components of Ha in the 
SN ejecta (Stathakis & Sadler 1991; Turatto et al. 1993). 
In the case of SN 1988Z, the underlying broad Ha com- 
ponent was also asymmetric, although not as severely as 
that of SN 2006tf, and it also increased in prominence 
with time (up to about day 115) and then faded again. 
In principle, though, the underlying SN ejecta would be 
easier to see in SN 1988Z because the total SN luminosity 
was lower than in SN 2006tf. 

5.3. Narrow Ha and HP Emission Profiles 

Our spectra at the four epochs are obtained with differ- 
ent spectral resolution (see Table 2), and this may be the 
primary reason for apparent differences in the narrow Ha 
profile (Fig.[9|).^ Only one epoch (day 66) had sufficiently 
high dispersion to fully resolve the line profile. It shows a 

"^ However, we do not believe that the lower spectral resolution 
on day 64 fully accounts for the weaker narrow absorption compo- 
nents of Fe II and other species, as compared to SNe 1994W and 
2006gy in Figure[2l because the spectrum with much higher resolu- 
tion obtained 2 days later does not show strong narrow absorptin 
features either. 
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Fig. 9. — Same as Figure |6] but highlighting the narrow Ha 
component. The scale has been adjusted so that the broad wings 
of Ha overlap. The day 66 spectrum (green) was obtained with 
DEIMOS using a higher resolution than the other spectra; it is 
shown at the observed resolution (solid green) , as well as smoothed 
(dashed green) to match the resolution on days 41 (black) and 95 
(orange). At degraded resolution, the Ha profile on day 66 is nearly 
identical to that on days 41 and 95, suggesting little or no change 
in the strength of the narrow component or its narrow P Cygni 
absorption. The day 32 spectrum (blue) was obtained with even 
lower resolution. The two vertical dashes mark the velocity of the 
P Cygni trough at —72 km s~^ and its blue edge at —125 km s~^. 

strong narrow emission component with a width of ~150 
km s^^ in both Ha and H/3 at the point where it meets 
the underlying broader base. This is an underestimate of 
the emission hue's intrinsic fuU width near zero intensity 
(FWZI) because a narrow P Cygni absorption feature is 
also seen in both lines, which partially absorbs the blue 
side of the line. Blueshifted absorption is stronger in H/3 
than in Ha, which was true for the broad components as 
wen (Fig.ig. 

The pre-shock CSM speed is best seen in the 
blueshifted P Cygni absorption feature of H/3 on day 66 
(Fig. [TO]), which has a blue edge indicating a progenitor 
wind speed of 190 km s"-'^. This is, again, similar to the 
case of SN 2006gy, where the blue edge of the narrow 
P Cygni absorption indicated a wind speed of roughly 
230 km s~^ (Smith et al. 2007). As noted in that paper, 
such wind speeds are much too fast for red supergiants 
and much too slow for Wolf-Rayet (WR) stars, but they 
are typical of blue supergiants and especially LBVs. As 
with SN 2006gy, this is likely to be a critical clue to the 
nature of the star that exploded and its immediate pre- 
SN mass loss that gave rise to its dense CSM. 

Despite the differences in spectral resolution, we find 
no clear evidence that the narrow Ha profile is changing 
with time. When we degrade the spectral resolution of 
the day 66 spectrum (the dashed green profile in Fig. [9]) 
to match the resolution on days 41 and 95, the narrow 
Ha emission component has the same profile shape and 
the P Cygni absorption feature is no longer seen. This 
suggests that the same narrow P Cygni absorption may 
actually be present on days 41 and 95 as well, but is not 
resolved in those spectra. The first epoch on day 32 has 
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Fig. 10. — Narrow component of the H/3 profile on day 66 from 
the high-resolution DEIMOS spectrum. This gives the best esti- 
mate of the progenitor's pre-shock wind speed of 190 km s~^ from 
the blue edge of the P Cygni absorption component. The minimum 
of the blueshifted absorption is at —80 km s~^. 

even lower spectral resolution, so a similar comparison is 
less clear. There is some suggestion of persistent narrow 
P Cygni absorption in H/3 at all epochs (Fig. [5]). 

The strength of the narrow Ha emission component 
relative to the underlying broad component changes 
slowly but steadily. If we measure the strength of the 
narrow emission compared to the total line fiux includ- 
ing the broader emission (listed as "N/T" in Table 2), we 
find that the narrow emission contributes roughly 17%, 
11%, 8%, and then 3% of the total Ha emission on days 
32, 41, 66, and 95, respectively (Table 2). As late as 
day 445, this fraction remains at about 3-4% (see be- 
low). Recognition that the narrow component makes a 
relatively minor contribution to the total line luminosity 
will be important in §5.5, where we discuss the behavior 
of the integrated line luminosities and equivalent widths. 

5.4. Late-Time Ha Emission 

Our late-time spectrum taken on day 445 is shown in 
Figure [11] This is the first epoch at which we see ev- 
idence that the overall character of the spectrum has 
changed significantly. The continuum emission is very 
weak, which is the reason for the very large equivalent 
width in Table 2. Besides Ha, faint emission from other 
Balmer lines and the Ca ii near-infrared triplet at 8500 A 
can also be seen. The Ha/H/3 flux ratio has now climbed 
to >11.6, indicating that it is dominated by collisional 
excitation rather than recombination (Raymond et al. 
2007). 

The broad components out to 7500 km s^^ (Fig-[8]) are 
no longer seen. This suggests that the material respon- 
sible for these features had reached the reverse shock by 
day 445. The absence of absorption, in particular, indi- 
cates that the power source of their underlying contin- 
uum had faded by this time as well, as expected if that 
continuum source was the underlying SN ejecta. 

The spectrum is dominated by a strong Ha emission 
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Fig. 11. — The late-time spectrum of SN 2006tf, obtained on 2008 Mar. 1, day 445 after discovery. The inset shows the Ho profile on day 
445, scaled arbitrarily and compared to the day 41 (shaded blue) and day 95 (orange) spectra from Figure|6] Narrow H/3 and [O III] A5007 
are real and may signify an underlying H II region, but there are many noise spikes redward of the broad near-infrared Ca II emission. 



line with a similar intermediate width as before, as well as 
a lingering unresolved narrow component that makes up 
a comparable fraction (about 3-4%) of the total emission 
as in the day 95 spectrum. The narrow component is 
about a factor of 2 narrower than before, though, with 
a FWHM of -80 km s"\ and reaching ±80 km s"^ at 
its base. The narrow P Cygni absorption component is 
no longer present, even though the late-time spectrum 
has higher resolution than our day 66 spectrum. The 
luminosity of the narrow component is a factor of '^6 
fainter than on day 95, with a luminosity of only 2.5 x 
10'^^ ergs s^^. This very narrow Ha emission is likely to 
be at least partially contaminated by a background H ii 
region, since narrow components of H/3 and [O in] A5007 
are also seen fFig. [TT|) . 

The inset of Figure [Tl] shows a detail of the Ha line, 
compared to the early-time profiles on days 41 and 95. 
The broader component of the line has a similar width as 
before, and shows complex behavior. At some velocities, 
the line wings match the symmetric day 41 profile, while 
at other velocities the late-time line wings are similar to 
those on day 95. These are probably important clues to 
the optical depth and geometry as the SN evolves. 

One very interesting quality of the late-time Ha profile 
is the pronounced asymmetry at relatively low velocities 
of ±1000 km s^^. Within this range, the intermediate- 
width component is very asymmetric and blueshifted as 
compared to the symmetric day 41 profile, with a sharp 
edge on the blue side where it meets the line wing. 

An intriguing possibility is that the net blueshift of the 
line may be due to dust formation in the dense post-shock 
cooling shell, analogous to the well-established case of the 
pecufiar SN lb 2006jc (Smith et al. 2008a), and possibly 
also the SNc Iln 2005ip (R. Chevalier 2008; priv. comm.) 
and 1998S (Pozzo et al. 2004). The optical depth would 
be highest at the limbs of that shell, favoring the extinc- 
tion of material at low redshifted velocities and speeds 
near zero. If dust has formed in the post-shock region 
of SN 2006tf, we would predict excess emission in near- 
infrared bands, analogous to the case of SN 2006jc, and 
possibly also signs of increased extinction. In that case, 
the fact that dust formed much later in SN 2006tf (days 



200-400 instead of by day 50 in SN 2006jc) can be under- 
stood as a consequence of the fact that SN 2006jc faded 
much more rapidly than SN 2006tf, and was much less 
luminous overall. Dust formation could also mean that 
the late-time continuum we measure is an underestimate. 
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Fig. 12. — The intrinsic luminosity of the Ha emission 
line in SN 2006tf, showing the total flux (solid dots), the 
broad/intermediate- width component from the post-shock gas (as- 
terisks), and the narrow component from the unshocked CSM (di- 
amonds), from Table 2. The decline rate of the narrow component 
of Ha in SN 2006tf follows the decline rate of the continuum lumi- 
nosity (the model diffusion curve from Fig. |4] divided by 2900; solid 
curve). For comparison, we also show the total Ho luminosities of 
SN 1994W (dotted hne) from Chugai et al. (2004), and SN 1988Z 
(dashed line) from Stathakis & Sadler (1991) and Turatto et al. 
(1993), as well as SN 1997cy. [These SN 1997cy measurements 
were taken from Fig. 11 of Pastorello et al. (2002). A. Pastorello 
(2008, private comm.) informed us that they were made from spec- 
tra presented by Turatto et al. (2000).] 



5.5. Line Luminosities and Equivalent Widths 



SN 2006tf Gets the Bronze 



13 



Figures [12] and [13] show how the strength of Balmer 
emission changes with time. We measured the equivalent 
widths of the total Ha and H/3 emission at the epochs 
for which we have spectra, and used the red continuum 
flux inferred from the Hght curve to derive the total Ha 
line flux.^ We also measured the Ha/H/3 flux ratio and 
the fractional contribution of the narrow CSM emission 
components to the total Ha flux in our spectra. These 
are listed in Table 2. Uncertainties in the Ha flux are 
typically 5-8%, depending on the noise in the continuum. 

The luminosity of the narrow component declines 
slowly, roughly following the same decline rate as the con- 
tinuum luminosity inferred from the light curve (Fig. 1121 
shows the same diffusion model curve from Fig. [3] divided 
by a factor of 2900 for comparison) . Thus, the narrow Ha 
emission represents a constant 0.034% of the bolomet- 
ric luminosity, and its effective equivalent width would 
be constant if it were measured relative to the underly- 
ing continuum. This implies that the narrow emission 
from the pre-shock CSM is radiatively excited by the 
same continuum source that powers the bolometric lu- 
minosity, but not the same source that powers the broad 
post-shock Ha emission. If this is Case B recombination 
emission, the mass of H gas needed to create the Ha line 
can be expressed as Mna — 11.4 M© (Lna/'^e), where 
Lfja is expressed in Lq and rie is the characteristic elec- 
tron density in the CSM in cm~'^. For the narrow Ha 
component in SN 2006tf, the value of Lho « 5 x 10^ 
Lq at early times implies a very large mass of pre-shock 
CSM — even for very high CSM densities, the total mass 
required is ~6 Mq (ne/10^ cm""^)"^. If this emitting 
material resides within a radius of ~10^^ cm, then the 
average density is ~10® cm~'^, implying that the progen- 
itor's wind had a mass-loss rate of ~0.1 Mq yr^^ for 
many decades before the SN, and prior to the even more 
extreme mass loss in the few years before explosion. 

The broad/intermediate- width component (and hence, 
the total Ha line luminosity) exhibits very different be- 
havior from that seen in the narrow CSM lines fFig. [TS]). 
Instead of fading with the continuum luminosity, the Ha 
luminosity of the post-shock gas rises dramatically as 
the SN fades. An important clue to this behavior is 
that the nature of the intermediate-width Balmer emis- 
sion changes dramatically during this time (see Table 
2): at early times, the Ha/H/? ratio matches the value 
one expects for recombination emission, suggesting that 
the heating of the post-shock gas emitting Ha is actu- 
ally dominated by photoionization heating. This changes 
with time, however, as the ratio climbs to Ha/H/3 >10, 
indicative of pure collisional excitation. It can also be 
seen in the equivalent- width behavior (Fig. I13|) . where 
Ha continues to rise, while H/3 levels off. 

This suggests that as time proceeds, SN 2006tf tran- 
sitions from an optically thick regime where radiative 
heating dominates the emission, to a regime where colli- 
sional heating by ongoing CSM interaction dominates. In 
some sense, this will be true for any strongly interacting 
SN Iln, because the highest optical depths are usually 
encountered at early times. One clear manifestation of 
such behavior is an increasing line-to-continuum ratio, 

* However, the Ha luminosity on day 445 was measured directly 
from the flux-calibrated spectrum, although that flux was uncertain 
by roughly 50% because of possible light cirrus clouds. 
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Fig. 13. — The equivalent widths of Ha and H/3 emission 
with time, measured in our spectra of SN 2006tf (emission-line 
equivalent widths are positive and include narrow, broad, and 
intermediate- width components). Evolution of the Ho equivalent 
width is shown for several other representative SNe Iln discussed 
in the text. Except for SN 1998S (Leonard et al. 2000), most ob- 
jects do not have published Ha equivalent widths, so we derived 
them from R-band photometry and published Ha line fluxes from 
several sources cited in the text. 



measured as the equivalent width. Figure [13] compares 
the Ha equivalent width of SN 2G06tf to that of several 
other SNe. All SNe Iln shown here have Ha equiva- 
lent widths that increase with time as the densities drop 
and the underlying continuum fades. However, it is clear 
that even though the Ha luminosity of SN 2006tf is not 
too different from some of these other SNe (Fig. [12]), its 
Ha equivalent width is systematically less than all oth- 
ers during the main peak of its light curve, and it does 
not catch up until late times (Fig. [T3)) . This discrep- 
ancy would be even more stark if we accounted for the 
likely fact that the explosion date should be at least 30 d 
before the discovery date, shifting the SN 2006tf data 
to the right in Figure [13] This is a sign that the contin- 
uum luminosity contributed by photon diffusion from op- 
tically thick material is much stronger in SN 2GG6tf than 
in the fainter SNe Iln, and that this component persists 
much longer. At early times, the fraction of the bolomet- 
ric luminosity contributed by the Ha line, iHa/iBoi, is 
about 0.4%, rising to about 2.6% at late times. Corre- 
sponding numbers for the other SNe Iln in Figure [T3l are 
Lna/ Lc ~ 2% at early times and about the same as the 
of SN 2006tf at late times. 

The different behavior with time in SN 2006tf suggests 
that it would be difficult to fully account for its high 
luminosity with a model like the one applied to most 
other SNe Iln, where Ha and continuum emission cool 
the post-shock gas in quasi-steady-state. The difference 
in the case of SN 2006tf is that the extra continuum 
luminosity may be due to delayed photon diffusion from 
an opaque shocked shell, as discussed earlier. Without a 
two-component model like this, it is difficult to see why 
the continuum luminosity of SN 2006tf would be a factor 
of 4-10 higher than that of SNe 1994W and 1988Z even 
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Fig. 14. — Spectropolarimetry of SN 2006tf on day 64. The top panel shows the observed polarization (estimated using a rotated Stokes 
parameter) in black with error bars, binned to 50 observed A pixel"^ for clarity. Overplotted in blue (shaded gray in the printed edition) 
is the total-flux spectrum from that date for reference. The continuum polarization between the lines is about 0.9%, with decreases in 
polarization visible at the wavelengths of the Balmer emission lines and especially O I A7774. The bottom panel shows the observed position 
angle of the polarization. 



though its Ha luminosity is only 1.5-2 times larger. 

Eventually, as the radius increases and the density 
drops, SN 2006tf becomes dominated by direct radia- 
tion from ongoing optically thin CSM interaction, very 
much like the late-time behavior of SN 1988Z and other 
SNe Iln. However, the drop in the total Ha luminosity 
at late times suggests that the progenitor mass-loss rate 
of SN 2006tf was not constant, but began to rise as the 
time of core collapse approached (see §7). 

6. SPECTROPOLARIMETRY 

The spectropolarimetry of SN 2006tf on day 64 is plot- 
ted in Figure [TH The continuum is polarized at about 
the 0.9% level, with several depressions in the polariza- 
tion at the locations of the strong Balmer emission lines 
and O I A7774. For Ha and O i especially, the depolar- 
ization seems to be primarily associated with the broad 
blueshifted components discussed in §5.2 and Figure [S 
The differing polarization of the lines and continuum in- 
dicates that at least some of the observed polarization 
must be intrinsic to the supernova and not simply due 
to ISP. We integrated the observed polarization over the 
spectral range 5050-5950 A in the rest frame of the su- 
pernova to simulate a rest-frame T^-band observation and 
obtained a value of Py — 0.91 ± 0.03% for the continuum 
polarization at a position angle of 6' = 135?4 ± 0?8. 

The few SNe Iln studied polarimetrically in the past 
have hinted that the objects as a class exhibit high polar- 
izations and hence large asymmetries. The first evidence 
that this might be true came from broad-band imaging 
polarimetric observations of SN 1994Y about 245 d after 
discovery (Wang et al. 1996). The i?-band polarization 
for that object (dominated by Ha) differed by more than 
1.5% (defined as the difference in the Stokes parameters 
q and u added in quadrature) from that measured in the 
B and V bands, which were dominated by continuum 
emission. Wang et al. (1996) concluded that SN 1994Y 
had significant intrinsic polarization, but without spec- 



tropolarimetry the interpretation was unclear. 

Subsequently, spectropolarimetric sequences of two 
SNe Iln have appeared in the literature. SN 1998S was 
observed within a week of explosion (Leonard et al. 2000) 
and twice at later epochs (Wang et al. 2001). Hoffman et 
al. (2008) recently presented three epochs of spectropo- 
larimetry of SN 1997eg ranging from 16 to 93 d after dis- 
covery. The initial continuum polarization was observed 
to be large in both objects (P « 2.0% in SN 1998S and 
P « 2.3% in SN 1997eg), with modulations larger than 
1.5% across the strongest emission lines. Both objects 
also showed time-variable polarization in the continuum 
and lines, with the polarization changing by more than 
a percent in the q — u plane. 

Despite the high data quality for these two SNe, in- 
terpretation has been complicated. Observed continuum 
polarization can be attributed to both ISP by dust along 
the line of sight and intrinsic continuum polarization of 
the SN due to electron scattering. Scattering by lines is 
believed to be intrinsically depolarizing, but if the line- 
scattering material is not distributed spherically, it may 
produce net polarization when the spatially unresolved 
SN is observed from afar. Finally, as discussed above, 
emission-line profiles show several components represent- 
ing material with potentially different spatial distribu- 
tions. The earliest epoch of data on SN 1998S (Leonard 
et al. 2000) showed separate polarization modulations in 
both narrow and broad lines, prohibiting any interpre- 
tation that allowed both components to be completely 
depolarizing. Hoffman et al. (2008) found that the po- 
larization of strong lines (Ha, H/3, and He i A5876) in 
SN 1997eg showed "loops" when plotted in the Stokes 
q — u plane. They attributed this to the differing geome- 
tries of narrow- and broad-line material, and thus to a 
mismatch in symmetry axes of SN ejecta and the CSM. 

Our observations of SN 2006tf do not have sufficient 
signal-to-noise ratio to disentangle all these effects, es- 
pecially from a single epoch of data, but we can make 
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TABLE 3 

Some Basic Physical Properties of SN 2006tf^ in a CSM Interaction Scenario 



Day 


LboI 


Tss 


Rbb 


-Rshcll 


c 


to 


McSM 


*Af 




(109 Lq) 


(K) 


(10^5 cm) 


(1015 cm) 




(10^** gcm-l) 


(Mq yr-i) 


(yr) 


32 


13.6 


7800 


4.52 


4.56 


0.98 


13.6 


4.1 


-7.6 


41 


12.7 


7500 


4.72 


4.72 


1.00 


12.7 


3.8 


-7.9 


64 


10.5 


6800 


5.22 


5.12 


1.04 


10.5 


3.1 


-8.6 


66 


10.3 


6800 


5.18 


5.16 


1.02 


10.3 


3.1 


-8.6 


95 


7.9 


6300 


5.26 


5.66 


0.86 


7.9 


2.3 


-9.5 


445 


0.64 


[6300] 


[1.5] 


11.7 


[0.016] 


0.64 


0.2 


-20 



a 



^ Kshell and properties in columns to its right arc derived assuming a constant shell expansion speed of 2000 km s 
in square brackets 

useful comparisons to SNe 1998S and 1997cg. The first 
point to note is that the observed continuum polariza- 
tion of SN 2006tf is lower than that of the other two 
objects, both of which showed polarizations of at least 
2% at early times. However, both objects showed signifi- 
cant ISP, making any estimate of the intrinsic continuum 
polarization difficult and model dependent. Also, the 
continuum polarization changed over time, with that of 
SN 1997eg decfining to ~1.5% by day 93. Our SN 2006tf 
observations were taken at least 64 d after explosion, so 
one possibility is that the polarization was higher at ear- 
lier times, but the optical depth to electron scattering 
declined with time, perhaps as a result of recombination 
or decreasing densities in the ejecta. 

Potential evidence that the intrinsic continuum polar- 
ization is actually lower in SN 2006tf than in the other 
two objects is that the depolarization of line features is 
significantly smaller. As mentioned above, the line fea- 
tures in the other objects were stronger than 1.5%. In 
SN 2006tf, the modulation at Ha is only about 0.4%. 
Quantitatively, the Stokes parameters {q,u) integrated 
over 6555—6575 A (the core of the narrow Ha line) are 
(-0.27, -0.60) with an uncertainty of ±0.09% for each. 
This value represents a difference of 0.41% in polariza- 
tion (Stokes parameters added in quadrature) from the 
y-band value quoted above, corresponding to ~l/4 of 
the polarization in the other objects. Overall, then, our 
general conclusion is that SN 2006tf is only moderately 
polarized at late times, showing no clear sign for extreme 
large-scale asymmetry in the CSM interaction region. 



Highly uncertain values are 




7. BASIC PHYSICAL PROPERTIES OF SN 2006tf 

In the preceding sections, we have described the basic 
energy demands, overall spectral properties, kinematics 
from line profiles, and line luminosities that are impor- 
tant clues to the nature and evolution of SN 2006tf. With 
these in hand, we can now quantify some of the funda- 
mental physical properties of the SN, listed in Table 3 for 
each day since discovery when we have spectra. Column 
2 lists the bolometric luminosity (Fig. [4]) derived from 
i?-band photometry, and column 3 lists the character- 
stic black-body temperatures for the continuum shapes 
in Figure [2l These are used to calculate the correspond- 
ing radius that a black body would have, listed as Rbb 
in column 4. In column 5, i?shcii is the radius of the shell 
expanding at a constant speed of 2000 km s~^, starting 
from a large initial radius suitable for Rbb at the early 
epochs. These two quantities are also plotted in Fig- 
ure [151 Column 6 lists the corresponding dilution factor, 
C = -^BsZ-'^sheii' which provides an estimate of the ef- 
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Fig. 15. — Filled dots show the derived black-body radius from 
Table 3, appropriate for the luminosity and temperature derived 
from the spectra in Figure [2] The solid line shows the radius for 
the observed nominal shock speed of 2000 km s"'^, and the two 
dashed lines show representative inner and outer radii where the 
thickness of the post-shock shell is ~15% of the radius. The dotted 
line shows radii for an expansion speed of 7500 km s~^ , where SN 
ejecta can almost reach the reverse shock by day 64 when the fast 
blueshifted emission is seen at those speeds. 

fective optical depth or geometric covering factor of the 
emitting regions of the shell. This dilution factor is a key 
component in the transition from optically thick to opti- 
cally thin, as we will discuss later. The next two columns 
give the wind-density parameter w — AfcsM/VcsM and 
the progenitor mass-loss rate Mcsm- These are the val- 
ues needed to account for the observed bolometric lumi- 
nosity of SN 2006tf in a simple CSM-interaction model 
with 100% efficiency, where the observed luminosity is 
emitted instantaneously by the post-shock gas, calcu- 
lated from the expression 



M, 



CSM 



= 2 L 



VcSM 



(1) 



where Vcsm = 190 km s~^ is the observed pre-shock 
CSM wind speed, and V^heii — 2000 km s~^ is the ob- 
served constant speed of the post-shock shell. Again, 
these provide an estimate of the mass-loss rate that cre- 
ated the CSM being swept up to produce the observed 
luminosity. Since this assumes 100% efficiency, it may 
be an underestimate; any clumping or global asphericity 
will raise the required values of w and Mcsm- Accounting 
for possible broadening of the line wings by electron scat- 
tering would raise Mcsm even further because it would 
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imply a smaller value of 14hcn- The last column lists the 
CSM flow timescale tj^j = i?sheii/VcsM, which is the num- 
ber of years prior to core collapse when the progenitor 
had the mass-loss rate in the previous column. 

The quantities in the last two columns are of partic- 
ular interest, because they provide a roadmap for what 
the progenitor star was doing in the decades before it 
died. This indicates that the progenitor underwent a 
sudden boost in mass loss in the decade just before ex- 
ploding. A caveat is that the instantaneous mass-loss 
rate may not have been quite as high as the peak val- 
ues around 4 M© , because diffusion of radiation from an 
optically thick shocked shell has the effect that the lumi- 
nosity from earlier CSM interaction, when the material 
was highly opaque, will effectively "pile up" and mimic 
a larger instantaneous mass-loss rate. Even allowing for 
this effect, though, average mass-loss rates of order 2 Mq 
yr~^ are needed in the decade before explosion. Before 
that time (at t minus 20 yr), the mass-loss rate was ap- 
parently 10 times less. 

As we have noted elsewhere for SN 2006gy (Smith et 
al. 2007, 2008b; Smith & McCray 2007), the only stars 
known to be capable of producing such extreme mass-loss 
rates of this order are LBVs during their giant eruptions 
(Smith & Owocki 2006), like the 19th-century eruption 
of rj Carinae when the star ejected 12.5 Mq in about a 
decade (Smith et al. 2003). Even the factor of 10 lower 
mass-loss rate of 0.1-0.2 M© yr~^ needed in the decades 
before that for SN 2006tf is still 1000 times stronger than 
the maximum rate that can be supplied by a line-driven 
wind of a massive star (Smith & Owocki 2006). How- 
ever, 0.1-0.2 M0 yr~^ is comparable to the less extreme 
LBV eruption of P Cygni in 1600 AD (Smith & Har- 
tigan 2006). Gal- Yam et al. (2007) have also discussed 
LBVs as potential SN Tin progenitors for independent 
reasons. If the stars that exploded as SNe 2006tf and 
2006gy were not actually LBVs, then they have done a 
remarkably good job of impersonating the H-rich com- 
position, amount of mass ejected, and ejection speeds of 
giant LBV eruptions. 

8. SUMMARY AND DISCUSSION 

8.1. Summary of Observational Clues 

The dataset we have presented above includes many 
different insights to the nature of SN 2006tf that weave a 
complex picture of a very luminous SN with remarkably 
strong CSM interaction. Figure [TBI shows an illustration 
of the essential structural elements of SN 2006tf follow- 
ing the hypothesis that it is a case of extremely strong 
CSM interaction, including diffusion from an opaque 
shocked shell, as mentioned many times throughout this 
paper. We now summarize the key observational clues of 
SN 2006tf presented in our study, along with the main 
implications of each and how they might be understood 
in a simplified model like that portrayed in Figure [TBI 

1. The total energy radiated at visual wavelengths is at 
least 7x 10^" ergs. It could be substantially more depend- 
ing on how soon after explosion SN 2006tf was discov- 
ered. If momentum is conserved, then at best, roughly 
half the initial kinetic energy can be converted into radia- 
tion in a CSM-interaction model. This suggests that the 
expanding shell should retain roughly an equal amount 
of kinetic energy. Indeed, our estimated shell mass of 
18 M© expanding at 2000 km s^^ does yield 7 x 10^" 



ergs. This requires a total initial explosion energy of at 
least 1.4 X 10^^ ergs (not including the kinetic energy 
of fast SN ejecta that have not yet reached the reverse 
shock). In other words, the very large CSM mass we es- 
timate is comparable to that which is needed to absorb 
momentum and decelerate the blast wave to only 2000 
km s^^. 

2. The decline rate of the optical light curve roughly 
matches that expected for the decay rate of ^^Co, in- 
cluding the late-time photometry more than 1 yr after 
explosion (Fig. |4]). The initial mass of ^^Ni needed to 
provide such a high luminosity is about 4.5 Mq. In that 
case, a pair-instability explosion would need to be in- 
voked, as in the case of SN 2006gy if it was also powered 
by radioactive decay (Smith et al. 2007). In our estima- 
tion, however, the strong signature of CSM interaction 
in the spectrum makes the radioactive-decay hypothesis 
less palatable for SN 2006tf. 

3. The decline rate of the optical light curve during 
its main peak is fit even better by a simple model of 
photon diffusion from an opaque shocked shell (Smith 
& McCray 2007). This model requires a shell mass of 
roughly 18 M© to account for the long diffusion time and 
high luminosity, which matches the large mass needed to 
decelerate the blast wave as noted in point (1) above. 
The initial radius of the opaque shocked shell is about 
180 AU or ~ 3 X 10^^ cm, close to the radius one might 
expect shortly before discovery (Fig.fTS]). 

4. However, both radioactivity and this diffusion 
model fall short of the luminosity needed to power the 
late-time photometry of SN 2006tf (Fig. S]). The similar 
CSM-interaction model by Woosley et al. (2007), requir- 
ing roughly 25 M© of CSM (but admittedly tailored for 
SN 2006gy instead of SN 2006tf), has the same difficulty 
(Fig. [5]). Thus, in any CSM-interaction model, we need 
to invoke additional, ongoing CSM interaction at large 
radii to explain the late-time behavior. Indeed, such on- 
going CSM interaction is evidenced by the strong Ha 
emission at late times. The progenitor mass-loss rate 
needed to account for the late-time luminosity through 
CSM interaction is about 0.2 M© yr~^, about 20 yr be- 
fore explosion. This is within a factor of 2 of the value 
required by the narrow CSM component (see next point). 

5. The narrow emission and P Cygni absorption com- 
ponents of Ha and H/3 reveal that the pre-shock CSM 
is expanding at roughly 190 km s~^. This is an impor- 
tant clue to the nature of the progenitor star, because 
it makes a red supergiant or WR-star progenitor seem 
unlikely, while the speed is comparable to those of LBVs 
and other blue supergiants. The luminosity of the nar- 
row Ha component implies a progenitor mass-loss rate at 
large radii outside the shock (^^10^^ cm; ejected about 20 
yr before core collapse) of at least 0.1 M© yr~^. This is 
close to the value needed to power the late-time luminos- 
ity (see previous point), and fully consistent within the 
uncertainty of the late-time luminosity estimate. This is 
also a factor of ~10 lower than the necessary mass- loss 
rate in the decade just before core collapse, signifying a 
sharp boost in Al immediately before the star's death. 

6. The intermediate- width component of the Ha line 
arises mostly in a swept-up, dense, post-shock cooling 
shell expanding at a constant speed of ~2000 km s~^ 
(Fig. [TB)) . This is the dominant speed of the forward 
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Fig. 16. — Cartoon illustration of the components of SN 2006tf at about 60 d after discovery, during the decline from the main light-curve 
peak. The primary feature is the massive post-shock shell of gas, composed of the swept-up opaque pre-SN envelope around the star 
ejected in the decade before core collapse. Most of the mass is in the cold dense shell (CDS), bounded by the forwrard shock (FS) and the 
reverse shock (RS). Diffusion of radiation from this shocked shell produces the main continuum photosphere (1) and the intermediate- width 
component of Hq. This shell expands at constant speed into the pre-shock CSM (dense wind of the progenitor). The interior of the shell is 
filled by freely expanding SN ejecta, the outermost parts of which are ionized by radiation (wavy lines) propagating inward from the reverse 
shock, exciting the broad He I and O I features seen in the spectrum. There is also a second photosphere (2) in the SN ejecta, which is 
fainter than the main photosphere and can only be seen if the main shell thins or develops clumps as time proceeds. The right panel shows 
a more detailed depiction of the post-shock gas, including the clumpy structure that forms due to instabilities in the cold dense shell layer. 
The dashed line here represents the photosphere at some arbitrary early time, working its way from left to right through the clumpy CDS 
as the SN expands. When it reaches a dense clump, the recombination photosphere will proceed through that clump, but for the regions 
between clumps it will eventually break through, allowing an observer to see the underlying SN ejecta. 



shock plowing into the CSM. This speed does not change 
perceptibly from day 32 onward. Since the shell does not 
decelerate even though it is emitting almost 10^^ ergs, the 
shell must already be very massive by day 32, consistent 
with our estimates above. 

7. The nature of the Balmer emission changes with 
time. At early times, the Ha/H/3 flux ratio is consistent 
with recombination, whereas at late times, the Hck/H/3 
ratio rises to more than 10, suggesting that it becomes 
dominated by direct coUisional excitation. 

8. Broad wings of Ha may be due in part to electron 
scattering, but there also appears to be an underlying 
broad emission component, seen almost exclusively at 
blueshifted speeds up to about —7500 km s~^ (Fig- 13). 
This broad component appears sometime after day 41, 
is seen on days 64 through 95, and disappears again at 
very late times. We propose that this feature corresponds 
to the outermost parts of the SN ejecta that have al- 
most reached the reverse shock (see Fig. [T6| . Material 
traveling at this speed would, in fact, just about reach 
the radius of the reverse shock by this time after explo- 
sion (Fig. [15]). The broad features are also seen in P 
Cygni absorption in He I A5876 and O i A7774. The ab- 
sorption requires some additional background continuum 
light source, which is likely to be the diffusion of radia- 
tion from the inner SN ejecta deposited by shock energy 



or radioactive decay. The luminosity required for the 
absorption strength implies that the underlying SN was 
overluminous as well, independent of CSM-interaction. 

9. A possible explanation for why the broad features 
are seen only from day 64 through 95 is that before that 
time, the shocked shell was highly opaque (the broad fea- 
tures reside interior to the reverse shock; Fig. I16p . Long 
after that time (by day 445), the SN ejecta luminosity has 
probably dropped far below that of the ongoing CSM- 
interaction region. 

10. The luminosity of the intermediate-width compo- 
nent of Ha is not correlated with the continuum lumi- 
nosity of the SN (Fig. [T^ . It rises as the continuum 
luminosity fades. Compared to other SNe Iln, the Ha 
equivalent width is lower, but rises to similar values at 
late times more than 1 yr after explosion (Fig. [T5)) . This 
is another clue that SN 2006tf has some additional source 
of continuum luminosity at early times, which is likely to 
be the slow diffusion of radiation from the massive swept- 
up opaque shell that mimics a normal H-recombination 
SN atmosphere, but at constant velocity. 

11. The intermediate- width post-shock Ha emission 
has pronounced asymmetry at late times, showing an 
asymmetric and blueshifted profile at velocities within 
roughly ±1000 km s^^. This may hint that dust forma- 
tion has occurred in the dense post-shock cooling shell. 
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analogous to SN 2006jc (Smith et al. 2008a), blocking 
gas mainly at the limbs and partially at the back side of 
the shell (bottom and top of the shell in Fig. \W\i . 

12. One epoch of spectropolarimetry of SN 2006tf re- 
veals relatively mild continuum polarization of 0.9% and 
also moderate depolarization of Balmer lines. Both ef- 
fects are weaker than in other SNe Iln for which data 
are available, so we conclude that SN 2006tf shows no 
severe, global asymmetry. Thus, the simplified spherical 
scenario we describe next is basically valid. 

8.2. Conclusion: Optically Thick to Optically Thin 

Altogether, then, we conclude that available evidence 
requires a two-component model for the power source of 
SN 2006tf. These two components correspond to differ- 
ent regimes of CSM interaction: (1) shock interaction 
with extremely dense and opaque CSM with progenitor 
mass- loss rates of order 1 M0 yr~^, where the mate- 
rial has such high optical depth that the radiation gen- 
erated by the shock interaction must diffuse out slowly 
through the massive amount of material after a delay 
corresponding to the diffusion timescale, and (2) more 
conventional, optically thin CSM interaction, where the 
luminosity generated by CSM interaction is able to es- 
cape without significant delay. As SN 2006tf expands and 
thins, it makes a transition between these two regimes. 

Exactly how and when does this transition occur? In 
a normal SN II-P, the H-recombination photosphere re- 
cedes backward through the expanding SN ejecta, main- 
taining a roughly constant Rbb during the plateau phase 
as it passes through ejecta with a large gradient in ex- 
pansion speed. The luminosity eventually drops sharply 
as the photosphere recedes all the way back through the 
H layer and the SN enters the nebular phase. This can- 
not be the case in SN 2006tf because the massive shell is 
expected to be geometrically thin and the observed ex- 
pansion velocity is constant, while we observe no sudden 
transition akin to SNe II-P. 

We propose that in SN 2006tf, the transition from an 
optically thick dense shell to the optically thin CSM in- 
teraction phase happened gradually as a result of clump- 
ing in the cold dense shell. The contact discontinuity, 
where material piles up between the forward and reverse 
shocks, will be severely Rayleigh- Taylor unstable and will 
not be a homogeneous constant-density spherical shell 
(Fig. [ini right panel). The size scale of density fluctua- 
tions will be similar to the thickness of the dense shell, 
which is geometrically very thin because it is relatively 
cold. Thus, the clumps that develop will have a typical 
size much smaller than the radius of the shell. The large 
number of small clumps across the surface of the shell 
will remain optically thick as recombination fronts move 
through them and their thermal energy is slowly radi- 
ated away, but the less dense regions between clumps will 
become optically thin more quickly, and the spacing be- 
tween clumps will grow with time until all the material 
in the shell eventually becomes optically thin. There- 
fore, we would expect that the observed transition from 
effectively optically thick to effectively optically thin can 
appear to be a smooth transition, where the shrinking ge- 
ometric covering factor of the optically thick clumps will 
mimic a shrinking average optical depth, while still al- 
lowing diffusion to be important for energy leakage from 
individual thick clumps. In effect, by looking through 



the increasingly porous shell, an observer can see the un- 
derlying SN ejecta even though the clumped material in 
the shell itself is still optically thick. 

This decrease in the effective covering factor of opaque 
chimps (rather than an actual decrease in the radius) is 
the "dilution factor" ( in Table 3, needed to explain the 
gradual fading as the shell expands. Thus, diffusion of 
shock-deposited energy from an opaque envelope domi- 
nates at early times (as in SN 2006gy; Smith & McCray 
2007), but as time progresses, the subsequent ongoing 
CSM interaction remains constant and eventually comes 
to dominate the late-time luminosity and spectrum as 
the optically thick shell fades. This is why the Ha/H/? 
flux ratio appears to rise steadily. 

At late times, this opaque shocked shell will have ra- 
diated away the bulk of its thermal energy and will cool 
significantly, heated only by the ongoing shock energy. 
The average densities in this 18 M© thin shell are ex- 
tremely high — roughly IQ-'^^ cm^"^ if the shell thickness 
is ^--^10% of its radius. With the low luminosity and large 
radius at late times, it would not be surprising if dust 
formation occurs in the opaque cool shell, for which we 
see suggestive evidence in the late-time blueshifted Ha 
line profile. We encourage near-infrared spectroscopy of 
SN 2006tf in order to determine if it shows the expected 
2.3 fim CO bandhead emission, as was seen in SN 1998S 
(Gerardy et al. 2000; Fassia et al. 2001). 

We might expect that this transition from thick to 
thin could, or should, happen for any SN Iln where the 
CSM close to the star is dense enough to be opaque, 
and this is probably true for other luminous SNe Iln at 
early times (e.g., SN 1994W; Chugai et al. 2004). The 
remarkable thing about SN 2006tf is just how high that 
luminosity was, how dense the CSM must have been, and 
how long the optically thick phase persisted. The wind 
density parameter, w, for SN 2006tf is orders of mag- 
nitude higher than that of most other SNe Iln. This is 
attributable to the very high mass- loss rate of the progen- 
itor of SN 2006tf, but also to the fact that its wind speed 
of 190 km s~^ is relatively slow (SN 1994W, for example, 
apparently had a CSM speed around 800-1000 km s~^; 
Chugai et al. 2004). One can therefore also understand 
the long duration of SN 2006tf in the context of its CSM 
speed: the terminal speed of 2000 km s~^ reached by the 
post-shock shell of SN 2006tf is still much greater than 
the 190 km s"^ speed of its CSM. In other SNe Iln with 
faster CSM, one might expect the interaction luminos- 
ity to drop as the shock decelerates, because the smaller 
difference in speed yields a weaker shock. 

The mass-loss rates of 1-2 M© yr~^ for about 10 yr 
before core collapse, and 0.1-0.2 Mq yr~^ for several 
decades before that, are truly astounding. The only stars 
known to be capable of producing it are very massive 
LBVs during their so-called "giant eruptions" (Smith 
& Owocki 2006). Normal stellar winds, even those of 
very massive and luminous stars, do not even come close. 
Thus, SN 2006tf adds a well established and rather ex- 
treme case to the growing body of evidence requiring 
massive eruptions in the decades preceding the core col- 
lapse of an initially very massive star. For SN 2006tf, 
the kinetic energy of that ejection was about 10^^'* ergs 
(assuming 18 Mq moving at 190 km s~^). In both SN 
progenitors and LBVs, these eruptions have no accepted 
explanation. For the SN progenitors, one potential ori- 
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gin conies from the pulsational pair instabihty discussed 
by Woosley et al. (2007) for SN 2006gy. That hypothe- 
sis, though, requires that the initial masses of these stars 
be 100-150 Mq, so any scenario would require that both 
SNe 2006gy and 2006tf must have marked the deaths of 
very massive stars. 

The type of diffusion model we suggest provides 
an explanation for one extremely important aspect of 
SN 2006tf where a traditional CSM-interaction model 
fails: the blast-wave speed of SN 2006tf remains steady 
at a surprisingly slow value of 2000 km s~^, even from 
very early times only 32 d after discovery, despite the fact 
that it is radiating almost 10^^ ergs in visual light. Tra- 
ditional CSM interaction predicts an observable deceler- 
ation with time (e.g., Chugai & Danziger 1994, 2003), 
while SNe 1988Z and 1998S provide observed exam- 
ples where deceleration is clearly documented (Chugai & 
Danziger 1994; Leonard et al. 2000; Pozzo et al. 2004). 
In the opaque shell-shocked model that we suggest for 
SN 2006tf, however, the bulk of the deceleration occurs 
very early as the blast wave runs through the opaque 



shell. Delayed photon diffusion allows the energy to be 
radiated away later. The massive circumstellar envelope 
that has been accelerated then coasts at a constant speed 
of 2000 km s^^ (so R/Rq is not large and adiabatic losses 
are minimized) while its thermal energy is radiated away 
as the shell expands, thins, and cools. 
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